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Abstract 

■ The link between the scahng function extracted from the analysis of (e, e') cross sections and the 
spectral function/momentum distribution in nuclei is revisited. Several descriptions of the spectral 

-T-^ function based on the independent particle model are employed, together with the inclusion of 

^ ' nucleon correlations, and effects of the energy dependence arising from the width of the hole states 

are investigated. Although some of these approaches provide rough overall agreement with data, 
they are not found to be capable of reproducing one of the distinctive features of the experimental 
scaling function, namely its asymmetry. However, the addition of final-state interactions, incor- 
porated in the present study using either relativistic mean field theory or via a complex optical 
potential, does lead to asymmetric scaling functions in accordance with data. The present analysis 
seems to indicate that final-state interactions constitute an essential ingredient and are required 
to provide a proper description of the experimental scaling function. 
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I. INTRODUCTION 



Investigations of inclusive quasielastic (QE) electron-nucleus scattering make it possible 
to obtain information about one of the main characteristics of nuclear structure, namely, 
the spectral function S{p, S) and its integral, the nucleon momentum distribution n{p) in 
nuclei l|-[3| . This provides insights into the validity of the mean- field approximation (MFA) 
and the role of the nucleon-nucleon (NN) correlations, as well as into the effects of Final-State 
Interactions (FSI) for inclusive electroweak processes. Using the shell model, it is possible 
in principle to obtain the contributions of different shells to S{p, S) and the momentum 
distribution for each single-particle state. However, due to the residual interactions, the 
hole states are not eigenstates of the residual nucleus but are mixtures of several single- 
particle states. This leads to the spreading of the shell structure and only a careful study 
of the momentum dependence of S{p, £) can separate the contributions from different shells 
(see, e.g. j4|). Such analyses have been carried out for few-body systems, complex nuclei 
and nuclear matter, focused mainly on the existence of high-momentum components of the 
nucleon momentum distribution due to NN correlation effects {2 7|. Since it is impossible 
within the MFA to describe simultaneously the density and momentum distributions in 



nuclei 
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12| . a consistent analysis of the role of the NN correlations is required using 



theoretical methods that go beyond the MFA to obtain a successful description of the relevant 



experiments. The present study uses the results found from studies of y-scaling 
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15| ) and superscaling (based on ■^/'-scaling variable, e.g. 15 



23|) obtained from analyses 
of inclusive electron scattering data. The latter consists in constructing a "superscaling 
function" f{il)) obtained by removing the single- nucleon content from the double differential 
cross section and plotting it versus a scaling variable '4){q, u). Scaling of the first kind of the 
scaling function {i.e., no explicit g-dependence of f{ip)) can be seen at excitation energies 
below the QE peak. Scaling of second kind [i.e., no dependence of f{'4') on the mass number) 
turns out to be excellent in the same region. When scaling of both first and second t ype s 



occur, one says that superscaling takes place. It was pointed out (see, e.g. [18 



20 



that the physical reason for the superscaling is the specific high-momentum tail of n{p ) 
which arises due to NN correlations and is similar for all nuclei. As was pointed out in 24| . 
however, a direct connection between the scaling function extracted from the analysis of the 
cross section data, and the spectral function only exists when one makes very restrictive 
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approximations. Along this line, caution should be kept in mind for the conclusions reached 
about the momentum distribution, because a close relationship between the latter and the 
scaling function also only emerges after some approximations are made. In particular, these 
are linked to the integration limits involved and the behavior of the spectral function 



In 



24j the analysis applied in the past to the scaling region (that is, negative values of the 



scaling variable y) was extended to positive y, leading to results that differ from those based 
solely on the scaling region and providing new insights into the issue of how the energy and 
momentum are distributed in the spectral function. 

Under certain approximations {e.g., see 2J] and references therein), in the case of plane 
waves the (e, e'N) differential cross section factorizes in the form: 

H PWIA 

= Ka'''{q,co;p,£,MS{p,£), (1) 



da 



de'dQ'dpj\[dQ 



N 



(e,e'N) 



where cx'^^ is the electron-nucleon cross section for a moving off-shell nucleon, S{p, S) is 
the spectral function that gives the probability to find a nucleon of certain momentum and 

In Eq. (dD p 



energy in the nucleus (see e.g. 



25H27j]) and is a kinematical factor 



28 



is the missing momentum and £ is the excitation energy that is essentially the missing 
energy minus the separation energy. Further assumptions are necessary 24 to show how 
the scaling function F{q,u) emerges from the Plane- Wave Impulse Approximation (PWIA), 
namely, the spectral function is assumed to be isospin independent and cx^^ is assumed to 
have a very mild dependence on p and £. The scaling function can be expressed in terms of 
the differential cross section for inclusive QE (e, e') processes: 

[da/de'dn%^^,^ 



Fiq,u) 



a'''^{q,u;p = \y\,£ = 0) 



(2) 



where represents the azimuthal angle-averaged single-nucleon cross section that also 
incorporates the kinematical factor K: 

A 



-eTV 



i=l 



a 



eNi 



2% 



Note that in Eq. a"^^ is taken a.t p = \y\, where the magnitude of the scaling variable y 
is the smallest value of the missing momentum p that can occur in the process of electron- 
nucleus scattering for the smallest possible value of the excitation energy {£ = 0), i.e., at 
the smallest value of the missing energy. Accordingly, in the PWIA the scaling function 



u) from Eq. ([2j) may be expressed in terms of the spectral function: 

F{q,u) = 2tx!! pdpdS S{p,S) , (3) 

J JT,{q,uj) 



where represents the kinematically allowed region (for details, see e.g. 2^). Only in 

the case when it is possible to extend the region S(g, w) to infinity in the excitation energy 
plane {i.e., at Smax oo), would the scaling function be directly linked to the momentum 
distribution of the nuclear system: 



oo 



n{p)= / d£S{p,S). (4) 







It was shown from the analyses of the inclusive electron-nucleus scattering that at high 
values of the momentum transfer the extracted scaling function Fc^p{q, u) becomes a function 



only of the scaling variable y, and not of q 



l|,B|l9 



19| . It was emphasized in |2J] that Eq. ([3 



does not apply to -Fcxp(5', oo) because of ingredients not included in the PWIA, such as final- 
state interactions, meson-exchange currents (MEG), rescattering processes, etc. 

Using the Relativistic Fermi Gas model (RFG) as a guide, the separate analysis of lon- 
gitudinal (L) and transverse (T) (e, e') data made it possible to introduce three "universal" 
experimental dimensionless superscaling functions: 

/exp(g, = A;Ai^cxp(g, uj); fc^Pig, uj) = kAF^J-piq, uj) , (5) 

kA being a phenomenological characteristic momentum scale for the specific nucleus being 
studied; this is equal to the Fermi momentum kp in the case of the RFG. In the present 
work we consider only the longitudinal scaling function and henceforth drop the subscript 
"L" for simplicity. Note that the effects of FSI and relativity on this function are often 



important and, as emphasized in |2J], any conclusion about the momentum distribution 
based on Eq. (|3]) should be made with caution. 

In the present work we study in more detail the relationship between the spectral function 
S{p,S) and the scaling function F{q,y). Our aim is to extract more information about 
the spectral function from the experimentally known scaling function, keeping in mind the 
restrictions of the PWIA. We take into account the effects of FSI and some other peculiarities 
of electron-nucleus scattering. We make an attempt to construct a spectral function that 
corresponds to the experimentally established scaling function following a series of steps 
on increasing complexity. Firstly, we construct S{p,S) within and beyond the Independent 



Particle Shell Model (IPSM). Secondly, we take into account FSI by computing the inclusive 
electron-nucleus cross section using the Dirac optical potential. We incorporate these results 
in the determination of the spectral function and consequently the superscaling function. 
In all steps we relate the results obtained for the scaling function to the empirical one. We 
establish a relationship between the single-particle widths obtained and the experimental 
ones. 

The theoretical scheme with a detailed analysis of the various approaches considered in 
the evaluation of the spectral function is presented in Sect. [Tll Here we also show and discuss 
the results obtained. We present a systematic study of the scaling function as well as the 
momentum distribution considering several different models and taking into account the role 
played by FSI. A summary of the work and our conclusions are presented in Sect. Illli 



II. SCALING FUNCTION IN RELATION TO THE SPECTRAL FUNCTION 
AND MOMENTUM DISTRIBUTION 



In this section we give the main relationships used in our approach in order to find 
a simultaneous description of the spectral function, momentum distribution and scaling 
function. As mentioned in the Introduction, the scaling function is given as a ratio between 
the inclusive electron-nucleus inclusive cross section and the electron-nucleon cross section 
at p = \y\ and £ = Within PWIA the scaling function is expressed in terms of the spectral 
function by Eq. ([3]). It was shown in 2J] that in this scheme the equations that relate the 
scaling function -F(g, y) with the spectral function in the regions of negative and positive 
values of the scaling variable y have the form: 
1 



y{q,y) 



ZTl 



d£S{p, £) 



if y < (6) 



r£~(p;Q,y) ry{q,y) 
pdp d£S{p, £)+ / pdp 

Je+{p;q,y) Jy 



£ {p;q,y) 



d£S{p, £) if y > . (7) 



In Eqs. (ED and ©: 



Y{q,uj) 



I + UJ)^J - M^^W - gAj /W^ 
A2-M0'iy2 + ^^|/yy2 



£^{p-q,uj) = {Ml + u) 



{q±pf 



m 



N 



M^'+p^ 



(8) 
(9) 



(10) 
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where u is the energy transfer, is the target nuclear mass, is the nucleon mass, 

ml + W^)/2 with 

^ in the form: 

^= 'r^ > ' (11) 



is the ground-state mass of the residual nucleus and A = {M^ 
W = \J (M^ + ujY — (f being the final-state invariant mass. 



In the RFG model the dimensionless scaling variable ip is introduced 15 

1 A-r 



^1 



y/(l + A)r + Kv/r(l + r) 



where rjp = kp/mN, = ^Jl + rjp — 1 is the dimensionless Fermi kinetic energy, k = 
g/(2mAr), A = uj/{2miy), and r = |(5^|/(4m^) = — A^ is the dimensionless absolute value 
of the squared 4-momentum transfer. The physical meaning of ip"^ is the smallest kinetic 
energy (in units of the Fermi energy) that one of the nucleons responding to an external 
probe can have. 



The scaling variables y and ip are closely related 



13, 



18| 



y_ 

kp 




m%l f y 



y_ 

kp 



(12) 



where vp is small, typically The dimensionless scaling function f{ip) is introduced 

{e.g. (isj) in the RFG model: 



/rfg(^) = kpF^Fci^) = - (1 - e (1 - ^2) . 



(13) 



As observed, the RFG model leads to a universal scaling function Jrpg which depends only 
on the scaling variable ip, but does not depend on the momentum transferred or on the 
nuclear species, that is, it superscales. 



A. Theoretical Spectral Functions: Independent Particle Shell Model and Beyond 

As noted in the Introduction the aim of the present work is to construct a realistic spectral 
function that leads to good agreement with the scaling function obtained from the inclusive 
electron-nucleus scattering data. We start with a given form for the spectral function, viz. 
that of the IPSM: 

SiPSAiip, ^) = E 2(2J. + iK{p)S{S - (14) 

i 

where ni{p) is the momentum distribution of the shell- model single-particle state i and 
Si is the eigenvalue of the energy of the state i. One may reasonably expect that when 
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effects beyond mean field are considered, the dependence on the energy here would be better 
represented by a function with a finite width in energy instead of by a (5-function. To explore 
the possible effect of a finite energy spread in Eq. (1141) . in what follows the energy dependence 
in Eq. f[T^ (the 5-function) is usually replaced by a Gaussian distribution G^X^ ~ ^d'- 

Sip,£) = Y,2{23, + l)nXp)G^X£-^i). (15) 

i 

where 

G,XE-E.>) = ^^e^ (16) 



and (Tj is a parameter for a given single-particle state i that is related to the width of the 
hole state i. In the present work we focus on the case of ^^O, together with a few results for 
^^C. For both nuclei we consider two parameters a\s and that are related to the widths 
of the Is and Ip hole states, respectively. We note that for simplicity we do not consider 
the differences in the spin-orbit partners 1^3/2 and lpi/2 states. We look for a best fit of 
the parameters in order to provide a good simultaneous description of the experimental 
scaling function /(V^), the experimental values for the widths of the hole states and the 
high-momentum tail of the momentum distribution. 

With a methodical purpose in mind we also consider and use for a comparison with 
the Gaussian in Eq. ( |T6i) another form of the energy dependence, namely the so-called 
"Lorentzian function" L-pX^ ~ ^i)'- 

1 Ti/2 

^^^^^~^'^ = ^{E-E,y + {VJ2Y' (17) 
where Fj is the width for a given single-particle hole state i. 

We start by taking riiij)) to be the momentum distribution of the harmonic-oscillator 
shell-model single-particle state i. Next, as can be seen from Eqs. ( fT5l) and ( fT6|) we account 
for the effects of nucleon correlations that give widths to the energy distributions of the hole 
strengths seen in (e, e') or (e, e'p) reactions. These widths may not in fact be symmetric, 
although in the present work for simplicity we limit ourselves to symmetric ones. In Fig. [1] 
are given the results for the scaling function compared with the longitudinal experimental 
data. There we have assumed equal values of and aip (cris = crip = cr) and vary cr in the 
region a = 10 — 90 MeV. In contrast, results in Fig. [2] and Fig. [3] correspond to fixed values 
of {Vis in the case of Lorentzian functions, see Eq. (fT7|) ) and aip {Tip) taken only at the 
extreme values 10 MeV and 90 MeV, and vice versa. 




FIG. 1: (Color online) Results for the scaling function f{ijj) for ^^0 obtained using HO single- 
particle wave functions (for a = 10 — 90 MeV) are compared with the longitudinal experimental 
data. The value of q is fixed to g' = 1 GeV/c. 

One can see that at fixed values of the parameter ais (Xis) and for running values of 
(Tip = 10 MeV {Tip = 10 MeV) and 90 MeV the main effects in the scaling function are 
observed in the shape of the curve and its maximum value which increases significantly as aip 
(Tip) goes down. Likewise the curve is extended in the negative ip region. On the other hand, 
in the case of fixed values of aip [Tip) and running values of ais = 10 MeV (Fi^ = 10 MeV) 
and 90 MeV, a smaller decrease of the maximum is observed and the discrepancies in the 
extended tail at < tend to disappear, giving rise to a similar shape for all ais (Tis) 
values considered. 

The main conclusion from our results presented in Figs. [TH3] is that, at least with a 
symmetric energy spread for the single-particle energy levels, it is not possible to get an 
asymmetry of the longitudinal scaling function similar to that shown by the data. 

Making use of Eq. (jl]) we calculate the IPSM momentum distribution which, as we are 
using HO single-particle wave functions, does not present a high-momentum tail. Our next 
step is to use natural orbitals (NOs) for the single-particle wave functions and occupation 
numbers employing a method where short-range NN correlations are taken into account. 
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a-is = 10 MeV (r^s = 10 MeV) is fixed = 90 IVieV (r^s = 90 IVieV) is fixed 




FIG. 2: (Color online). Results for the scaling function /{t/j) for ^^O obtained using HO single- 
particle wave functions. The values of ais {Tig in the case of Lorentzian functions, see Eq. ()17p ) 
are fixed and aip = 10 MeV (Tip = 10 MeV) and dip = 90 MeV {Tip = 90 MeV) have been used. 
The results are compared with the longitudinal experimental data. The value of q is fixed to q = 1 
GeV/c. 



In what follows we use the NO representation of the one-body density matrix (OBDM) 
obtained within the lowest-order approximation of the Jastrow correlation method 291 ]. 

The NOs v^a(r) are defined 30| as the complete orthonormal set of single-particle wave 
functions that diagonalize the OBDM: 

p(r,r') = 5^iV,v.:(r)v..(r'), (18) 

a 

where the eigenvalues Na (0 ^ Na ^ 1, J2a = A.) are the natural occupation numbers. 

The NO single-particle wave functions are used to obtain the momentum distributions 
ni{p), and from them the spectral function according to Eq. ( fT5i) . The results for the scaling 
function obtained using NOs and HO single-particle wave functions for various values of the 
parameters ais and dip are given in Fig. HI They are represented by solid (HO) and dashed 
(NO) lines and compared with the RFG result (dotted) presented also for reference. As 
observed, the main effect introduced by the use of NO, compared with HO, is an enhancement 
in the maximum of the order of ~10%. In contrast, the tail is slightly reduced in the region 
of negative ip. However, notice that both models lead to a weak asymmetry in the scaling 
function f{ip) that is not in accordance with the significant tail extended to positive ip- 
values, seen in the analysis of (e, e') data. It should be mentioned here that in addition to 



a-ip = 10 MeV (r^p = 10 MeV) is fixed 



Oip = 90 IVIeV (r^p = 90 IVIeV) is fixed 
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FIG. 3: (Color online). Results for the scaling function /{i/j) for ^^O obtained using HO single- 
particle wave functions. The values of dip (Tip) are fixed and ais = 10 MeV (Fig = 10 MeV) and 
ais = 90 MeV (Fi^ = 90 MeV) have been used. The results are compared with the longitudinal 
experimental data. The value of q is fixed to g = 1 GeV/c. 



Is- and Ip- components there are also Id- and 1 /-components in the NOs obtained from the 
Jastrow correlation method. Unless specified otherwise we take the same values of a for all 
of them. 

In Fig. |5] we present the evolution of the scaling function /{ip) for different values of 
q running from 100 to 2000 MeV/c. Results have been obtained making use of the HO 
momentum distributions for the Ip- and Is-shells in ^^O. From these one gets the spectral 
function according to Eqs. (fTSl [T6l) and finally the scaling function using the expressions in 
Eqs. (El [7j). As already mentioned, the HO model is not capable of producing the strong 
asymmetry observed in the data, but it can be seen that for q > 600 — 700 MeV/c scaling 
of first kind is fulfilled. 

In what follows we focus on the study of the momentum distribution obtained in different 
approaches. Results are presented in Fig. [6]where we use a log scale in order to emphasize the 
differences. The RFG momentum distribution is compared with the result obtained using 
HO single-particle wave functions and the NOs from the Jastrow correlated approach 29 |. 
The nucleon momentum distribution n^^^ obtained in 23|, ISlJ by using the Light-Front 



Dynamics method (LFD) 32|is presented in Fig. [6] as well. We also show in the same 
figure the results obtained in 33|| within the Relativistic Mean Field (RMF) model with and 
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-2-1012 

FIG. 4: (Color online) Results for the scaling function f{ip) for ^^O obtained using NOs (dashed 
lines) and harmonic oscillator (HO) single-particle wave functions (solid lines) for various values 
of the parameters ais = crip = 10 MeV (red lines), 30 MeV (green lines) and 70 MeV (black lines). 
The RFG results are shown for comparison (dotted line). The value of q is fixed to g' = 1 GeV/c. 



without taking into account FSI. Finally, we give in Fig. |6]as an example (being considered 
and used in our previous work) the result obtained with the Coherent Density Fluctuation 
Model (CDFM) (Q, S, 11|; see also 21, 23|). The model is a natural extension of the RFG 
model. It is based on the 5-function limit of the generator coordinate method {34 1 and 
accounts for long-range NN correlations of collective type. In it the nucleon momentum 



distribution has the form: 



where the CDFM Wigner distribution function is: 



n(k) = / drW{r,k), 



(19) 



W{r,k) 



with 



and 



Vr,.(r,k) 



kpix) 



(27r)3 
37r2 



Q{x - \v\)e{kF{x) - |k|) 



1/3 



3/1 

Atix'^ 



(20) 



(21) 



(22) 
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FIG. 5: (Color online) The scaling function /{i/j) for ^^O calculated using HO single-particle wave 
functions (with o" = 0) is presented as a function of ip for a range of fixed values of q extending 
from 100 MeV/c (right-most curve) to 2000 MeV/c (left-most curve) in steps of 200 MeV/c. For 
q > 700 MeV/c scaling of the first kind is seen to be fulfilled. 



The model weight function F{x) is obtained by means of a known density distribution p{r) 
for a given nucleus: 

1 dp{r) 



\F(x) 



dr 



(at dp{r)/dr < 0). 



(23) 



The CDFM has been applied to studies of the superscaling phenomenon (e.g. 20N23|) in 
inclusive electron-nucleus scattering, as well as to analyses of neutrino and antineutrino 
scattering by nuclei of both charge- changing 23|, l35|] and neutral-current 36|] types. 

Before entering into a discussion of the RMF results, it is interesting to point out the 
significant differences introduced by the other models. As noticed, the presence of nucleon 
correlations (through the Jastrow approach) leads to a significant tail at high momentum 
values compared with the pure HO result. On the contrary, the momentum distribution is 
slightly reduced for small-to-intermediate p. As can be seen, the tail at high p is extremely 
large for n^™(p) and is also present in the case of the CDFM model. The latter is a result 
of the NN correlation effects accounted for in the model. Later in this work we also show 



the role of these correlations within the CDFM approach 23|] (with phenomenologically 
introduced FSI effects) on the scaling function f{ip) in the cases of ^^C (Fig. [9]) and ^^O 
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(Figs. [H and [n]). 
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FIG. 6: (Color online) Results for the momentum distribution for ^^C obtained using harmonic 
oscillator single-particle wave functions (HO - green line); NO taken from the Jastrow model 
(Jastrow - blue line); Relativistic Fermi Ga s ( RFG - red line); CDFM results (CDFM - grey 



area 
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23|); LFD results (LFD - olive line [23|). Normalization: / dkn(k) = 1. 



Next we discuss in more detail how the momentum distribution function is evaluated 
within the general framework of the RMF. The particular case where FSI are turned on is 
considered in the next subsection. In the RMF case without FSI the momentum distribution 
fiRMF^p^ is given in the form: 

n^MF^p^ = J2 ^«heu(p) (24) 

shells 

with 

nshell(p) = I '/'shell (P) I ^ (25) 

where v^sheii(p) are single-particle 4-component wave functions that are solutions of the Dirac 
equation with an RMF (relativistic Hartree) potential. 

It is also interesting to consider the momentum distribution for the case in which the 
relativistic bound nucleon wave functions are projected over positive energy components. In 
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this case the four-spinors v^shci^p) go into four-spinors fshcwip) by the action of the positive- 
energy projection operator 

1 </?shcii(p) • (26) 
Hence the "positive-energy projected" momentum distribution results: 

n^MF^p^ = Yl = Yl l^sheiib)!' . (27) 

shells shells 

The latter from Eq. fl27|) is labeled in Fig. [H] as ''no FSI, projected" , although we should 
point out that it differs from the usual non-relativistic analyses where two-component 
Schrodinger-like equations are considered and a non-relativistic reduction of the relativistic 
4x4 current operator into a 2 x 2 form is used. Here, we project out the negative-energy 
components of the bound nucleon wave functions while a fully four-component description 
of the Dirac spinors is maintained. The outgoing nucleon is described by means of a rel- 
ativistic (Dirac) plane wave. These two assumptions, i.e., plane waves in the final state 
and positive-energy projection in the initial one, lead to the so-called Plane- Wave Impulse 
Approximation (PWIA), which should be distinguished from the Relativistic Plane- Wave 
Impulse Approximation (RPWIA) where the plane-wave approach for the ejected nucleon 
is also considered, but where the fully relativistic bound nucleon wave function (without 
projection) is used. 

In PWIA the exclusive electron (e, e'N) cross section factorizes in the form: 

'^'^ = Ka'^'ustenip) , (28) 

PWIA 



dn'de'dn 



N 



where K is a. kinematical factor and a'^^ is the electron-nucleon cross section. Thus Eq. ( 128|1 
makes it possible to obtain the momentum distribution for different shells by means of the 
cross section calculated within the PWIA: 



?^shcll(j5) 



da 




_dVL'de'dVLN_ 


PWIA 



(29) 



This result coincides with the one given through the positive-energy projected wave functions 
in Eq. (ETj). 

Results for n^^^{p) (labeled as "no FSI, RMF") and n^^^(p) (no FSI, projected) are 
compared in Fig. [61 As shown, n^^^(p) follows the HO result closely, showing a steep slope 
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with increasing p. Likewise, n^^^{p) presents a similar behavior but without the diffraction 
minimum due to the presence of the lower components in the relativistic bound nucleon 



wave functions. These results confirm previous studies presented in [37 1. 



B. The role of FSI 

The factorized form of the (e, e'A^) cross section shown in Eq. ( l28l) does not apply when 
final-state interactions are included in the description of the reaction mechanism. However, 
one can introduce a function given as the ratio between the differential cross section and the 
single-nucleon cross section (multiplied by the kinematical factor K). This is usually called 
the reduced cross section or distorted momentum distribution: 

da 



Pip) = n'^'^'ip) 



dQ'de'dn 



N 



FSI 



(30) 



Ka^^ 

Although in Eq. (!30|) we show the function n'^'*** to be only dependent on the missing 
momentum p, this is not so in general because of the non-factorized form of the (e, e'A^) 
differential cross section when FSI are included. The distorted momentum distribution de- 
pends not only on the missing momentum p but also on the other independent variables 
in the scattering process, for example the momentum transfer q. Hence, calculations for 
different kinematical situations may lead to slightly different distorted momentum distribu- 



tions. However, this dependence has been proven to be weak in most cases |38|-|40[| . The 
resulting effective momentum distribution n'^"^*(p) for a momentum transfer of 1 GeV/c and 
the momentum of the final nucleon equal to the momentum transfer is also displayed in 
Fig. [6] (labeled as "FSI, RMF"). In computing this effective momentum distribution using 
Eq. (!30|) . the same real scalar and vector potentials that describe the bound nucleons are 



4l| has 



employed to describe the final nucleon. The nonlinear effective interaction NLSH 
been used. 

As observed, FSI evaluated within the RMF approach gives rise to the presence of a very 
significant tail at high missing momentum which is several orders of magnitude larger than 
the results obtained within the plane- wave approach. On the contrary, n'^"^*(p) is smaller for 
p-values below the Fermi momentum. It is important to point out that the right amount 
of asymmetry in the longitudinal scaling function required by (e, e') data is also reproduced 
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by a semi-relativistic shell model when the continuum states are described with the Dirac- 



42, 43l- 



equation-based (DEB) equivalent potential to the RMF employed here 

Once the distorted RMF momentum distribution has been obtained [Eq. ( I30p ]. we cal- 
culate the "distorted" spectral function S'^^^^{p,S) using Eqs. f[T5]l and flTBl) . Finally, the 
corresponding scaling and superscaling functions can be evaluated making use of the proce- 
dure outlined through Eqs. ([5l [6l [7]). The results for f{ip) are given for different values of 
cr = ais = (Tip (from 10 to 100 MeV) in Fig. [71 Also, for completeness, we show in Fig. [7] 
the scaling function fXip) obtained within the RMF+FSI model, but evaluated directly from 
the inclusive (e, e') differential cross section [33|, that is, 

RMF+FSI 



da 



de'dn' 



(e,e') 



(31) 



''^ a''^{q,u]p= \y\,£ = 0) ' 

Results in Fig. [7] show that the RMF with FSI included gives rise to an asymmetrical scal- 
ing function which follows the behavior of data, irrespective of whether we use Eqs. fll5|30p 
or we use Eq. ( 13T]) . In the first case, significant discrepancies emerge when different values of 
the (T-parameter are used in the Gaussian energy-dependent functions. The scaling function 
obtained from Eq. fl3ip is close to that obtained from Eqs. (11511301) for cr = 30 MeV. We 
note that for the negative ip region, the scaling function obtained from the RMF calculation 
for the cross section lies well below the data. The comparison with data in this region is 
thus improved when a finite width for the energy is introduced in the spectral function. As 
observed, the case a = 30 MeV is the one that agrees better with experiment and with the 
theoretical curve obtained directly from the inclusive cross section (see Eq. (|3T]) ). We use 
a = cTis = (Tip to check the dependence of the scaling function f{ip) on cr. Below we dis- 
cuss the optimal values of these parameters ais ^ aip, and check the results for the scaling 
functions when these parameters give widths close to the experimental ones. 

In addition to the description of FSI within the framework of the RMF model, we have also 
considered alternative analyses of final-state interactions in the case of inclusive electron- 
nuclei processes. In what follows we discuss a method proposed in 4^ and present the 
results we obtained with this method. 

According to the authors of {4^] two types of FSI effects, Pauli blocking and the interaction 
of the struck nucleon with the spectator system described by means of the time-independent 
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optical potential (OP) 



U = V-iW 



(32) 



proposed in 



45|, can be accounted for by replacing in the PWIA expression for the inclusive 



electron-nucleus cross section 



dat 



the energy-conserving delta-function by 



EpEp/ 



6{u + M-E 



6{u + M-E-Ep,)L';:^H, 



W/tx 



em, t 



(33) 



(34) 



[io + M-E-Ep,-V]^' 
In Eq. ( I33ll the index t denotes the nucleon isospin, L^™ is the leptonic tensor, H^!^^ j is the 
hadronic tensor and St{p,E) is the proton (neutron) spectral function. The terms E]f_, Ep, 
Epi and E represent the initial electron energy, the energy of the nucleon inside the nucleus, 
the ejected nucleon energy and the removal energy, respectively (see 4J] for details). 

The real and imaginary parts of the OP in Eqs. ( I32]l34l) are obtained from the Dirac OP 
in 46|. We use spatially averaged values, evaluating them at the r- values that match their 




FIG. 7: (Color online) Results for the scaling function f{ip) obtained using single-particle wave 
functions taken from RMF + FSI for calculation of ni{p) and for various values of the parameters 
a = (Jis = aip. Also shown are results obtained through Eq. ()3ip . Momentum transfer fixed to 
q = 1 GeV/c. 
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FIG. 8: (Color online) The real and irnaginary parts of the optical potential U{p') for ^^C calculated 
by Eq. ([35]) using a Dirac OP from |4a] and the imaginary part of the optical potential for ^^O 
given by Eq. ([36]) . 



respective root mean squared radii 46)]. As a result the OP U{p') related to the scalar (S) 
and vector {V) part of the potential in 46|| is obtained in the form (see also 44l|): 



Ep/ 



U{p') = J[M + S{Tp,,fsW + P" + V{Tp,,fv). 



(35) 



We consider also an OP with the imaginary part of the potential U{p') given in |47j : 



W 



he 



IP 



—Pnucl<^NN-„ 



(36) 



where the values of pnuci and (Tattv for ^^O are taken to be: pnuci = 0.16 fm~^ and (Tattv = 
40 mb. 

In Fig. |H]we give the calculated real and imaginary parts of the OP f/(p') from Eq. ([53]) 
and the imaginary part of the OP from Eq. (1361) . As noted in 4^, for |p'| > 580 MeV/c 
the real part of t/(p'Iis positive, which is inconsistent with correlated Glauber theory 48 |. 



Therefore, following 
was noted in Ref. 



44j when |p'| > 580 MeV/c the real part of the OP is set to zero. As 
the real and the imaginary part of the optical potential quantitatively 
have very different effects, and describe different aspects of the FSI. The imaginary part 
W accounts for two-body scattering processes involving large momentum transfers, which 
lead to a strong damping of the motion of the recoiling particle, whereas the real part \J 
produces a shift of its energy. The effect of the imaginary part is known to be dominant at 
large momentum. 
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FIG. 9: (Color online) Results for the scaling function /(V') with and without taking into account 
FSI in the RFG model (for a given momentum transfer q = I GeV/c and energy of the initial 
electron e = 1 GeV) are compared with the longitudinal experimental data and CDFM results 



using ci =0.63 



23|. 



Following the approach of 



44l | we calculate the scaling function according to Eq. 



but evaluating the cross section in the numerator by accounting for FSI using the Dirac 
OP, i.e., exchanging the (5-function in Eq. (133|1 as shown in Eq. f lM|) . The expression for 
S{p,S) in Eq. (155]) is taken in the form in Eq. fllSp . where we use momentum distributions 
ni{p) obtained in different approaches. The results in the case of the RFG model spectral 
function are presented in Fig. [9] for a given momentum transfer q = 1 GeV/c and energy of 
the initial electron e = 1 GeV with and without accounting for FSI. The scaling function 
f{ip) using the Sipsm{p,^) of Eq. f|T4|) with HO single-particle wave functions is presented 
in Fig. [To] for q = 1 GeV/c and e = 1 GeV. In this case we consider two different types of 
the time-independent optical potential for ^^O: i) the one obtained by Eqs. fl32|34|35p using 
the scalar and vector part of the potential from 46| and ii) using the imaginary part of the 
potential U{p') from Eq. ([36]) [47 1. 

As can be seen from Figs. [9] and [10] FSI lead to a decrease in the maximum of the scaling 
function, while an enhancement of the tails is observed for both negative and positive values 
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FIG. 10: (Color online) Results for the scaling function f{ip) with and without taking into account 
FSI using the spectral function in the IPSM (Eq. (jl4p ) with HO single-particle wave functions 
(for a given momentum transfer q = 1 GeV/c and energy of the initial electron e = 1 GeV) are 



compared with the longitudinal experimental data and CDFM results using ci = 0.63 



2M- 



of tp. Also from Fig. [TO] can be seen the important role of the type of the time-independent 
optical potential used. 

In Figs. [9] and [101 as well as in Fig. [11] we give for a comparison the results for the 
scaling function obtained within the CDFM approach 23| in which FSI effects are intro- 
duced phenomenologically. As can be seen, accounting for FSI together with the effects of 
NN correlations is necessary for a reasonable explanation of the experimental data for the 
longitudinal scaling function, especially the maximum and the tail of f{ip) aX ip > 1. 

Finally we consider results that correspond to a spectral function constructed with NOs 
from the Jastrow correlation method, within this alternative way to take into account FSI. 
These results are presented in Fig. [TT] for q = 1 GeV/c and electron beam energy e = 
1 GeV (we consider two different time- independent optical potentials for ^^O, the same as 
in Fig. [To]) . The value of the parameter cti^ is fixed to be 8.7 MeV (that corresponds to the 
experimental width of the Is-state in ^^O taken from j49]) and aip = aid = ^i/ = 0-5 MeV. 
The results depend very weakly on the choice of parameters cTj. 
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FIG. 11: (Color online) Results for the scaling function f{ip) taking into account FSI using NOs 
from the Jastrow correlation method (for a given momentum transfer q = 1 GeV/c, energy of the 
i„i.ial ele*„„ . ^ . GeV, a„d pa..„.ete. ^ B.7 MeV, ^ ^ ^ 0.5 MeV) 
compared with the longitudinal experimental data and CDFM results using ci = 0.63 |23l ]. 



The next step of our work was to obtain optimal values of the parameters cTj in the 
Gaussians Gf^.{S — Si) [Eq. (fT6|) ] by an additional fitting procedure. These values of cTj are 
compared with the experimental widths of the hole states. The results of our calculations 
in the case of the RMF model accounting for the FSI's are given in Fig. [T21 Here the value 
of the parameter ais is fixed to be 7.8 MeV (that corresponds to the experimental width of 
the Is-state in ^^C taken from 49[) and the values of aip are varied. In Fig. [12] is given also 
the result for the scaling function obtained when more recent experimental data for the Full 
Width at Half Maximum (FWHM) for Is- and Ip-states (taken from 50|) are used. Notice 
that the latter result (black dotted line) almost coincides with the case uij, = 5 MeV and 
(Tis = 7.8 MeV. As can be seen, the results for the scaling function cannot reproduce very 
well the experimental data using experimental widths. 

Finally, in Fig. [13] are presented and compared the results for the scaling function f{ip) 
calculated using single-particle functions taken from HO and RMF+FSI and for two energy- 
dependence functions, Lorentzian and Gaussian. We note the good agreement of f{ip) (using 
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FIG. 12: (Color online) Results for the scaling function /(V') obtained using single-particle wave 
functions taken from RMF + FSI for calculation of ru(n). Five curves are shown. Four (given by 



solid lines) correspond to ais = 7.8 MeV taken from 



and aip ranging from 5 MeV (green line) 



to 20 MeV (blue line) in steps of 5 MeV. The fifth curve (black dotted line) corresponds to widths 



for Is- and Ip-states taken from 
transfer: q = 1 GeV/c. 



50|] (FWHMi^ = 20 MeV and FWHMip = 6 MeV). Momentum 



widths close to the experimental Tip = 6 MeV, Fi., = 20 MeV and single-particle functions 
taken from RMF+FSI), including the maximum value and the asymmetry of the scaling 
function. In this way, together with the results for the momentum distribution and the 
scaling function, we obtain a completion and self-consistency of our study. 



III. CONCLUSIONS 



The analysis of the world (e, e') data has clearly demonstrated the validity of scaling 
arguments. Not only has the scaling phenomenon been proven to be fulfilled, but also 
the specific shape of the scaling function, with a significant tail extending to high values 
of the scaling variable, has emerged from the separated (e, e') longitudinal data. These 
results constitute a strong constraint for any theoretical model describing electron scattering 
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FIG. 13: (Color online) Results for the scaling function /(V') obtained using single-particle wave 
functions taken from RMF + FSI for calculation of nj(p) (parameters Pip = 6 MeV, Fi^ = 20 MeV 



are fixed to the experimental widths of the Ip- and Is-states in ^^C 
functions, Lorentzian and Gaussian. 



soi l) and two energy-dependence 



reactions. Moreover, the scaling function is an observable which has been used in the past 
in order to get information on the spectral function S{p,£) and/or the nucleon momentum 
distributions of nucleons in nuclei n{p). However, as explored in detail in our previous 
work 2J], the connection between the scaling function (given directly from the analysis of 
(e, e') data) and S{p,S) or n{p) only exists under very restrictive conditions: i) the Plane- 
Wave Impulse Approximation (PWIA) in the description of the reaction mechanism, and ii) 
additional assumptions on the integration limits consistent with the kinematically allowed 
region. Being aware of these restrictions and with some caution in the discussion of results, 
the link of the spectral function/momentum distribution with scattering observables remains 
a delicate issue in nuclear physics. Effects beyond the PWIA, such as final-state interactions 
(FSI), meson-exchange currents, rescattering processes, etc., which may have very important 
roles to play in describing (e, e') processes, should be incorporated in the general analysis. 



In this paper we extend the previous work presented in 



241 ] by considering several models 



dealing with the nuclear structure, and including as well different descriptions to account 
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for FSI. From the momentum distributions provided by three nuclear models: Harmonic 
Oscillator, Natural Orbitals with Jastrow-type correlations and Relativistic Mean Field, 
we construct the spectral functions where the energy-dependence given through the delta 
function has been changed by introducing Gaussian distributions. In this way, we present 
a systematic analysis for different values of the Gaussian distribution parameters related 
to the width of the hole states. Finally, the scaling function is evaluated by making use 
of the general derivative expressions that connect it with the spectral function. We have 
explored in detail the role played by Gaussian or Lorentzian energy dependences in the 
scaling function showing that, although significant differences may emerge, all models based 
on PWIA lead to scaling functions which do not show the asymmetry that data exhibit, 
with an extended tail at high positive values of the scaling variable ip. At the same time we 
note that the width in energy introduced for the hole state, as opposed to the (5-function of 
the strict IPSM, can however improve the agreement with data in the negative ip region of 
the scaling function. 

Effects of FSI have been proven to be essential in the description of electron scattering re- 
actions. Hence, we have included FSI in our analysis within the framework of the RMF and 
the method using a complex optical potential {4^. In the former we obtain the reduced cross 
section (or effective distorted momentum distribution) for each shell as the ratio between 
the exclusive (e, e'A^) cross section evaluated within the RMF model, and FSI included, and 
the corresponding single-nucleon cross section. These "momentum distributions" incorpo- 
rate the role of FSI, and they can be used to get the spectral function following the same 
procedure as in the previous case (PWIA). Again, we perform a systematic analysis by con- 
sidering different values of Gaussian and Lorentzian widths. Finally, the scaling function is 
directly evaluated from the spectral function. We have compared these results with the one 
obtained directly from the ratio between the inclusive (e, e') cross section (evaluated within 
the RMF+FSI model) and the single-nucleon cross section taken at p = \y\ and £ = 0. In 
all cases the results obtained clearly show the essential role of FSI in producing the required 
asymmetry in the scaling function, i.e., to be consistent with data. Concerning the second 



method we explored for including FSI following |4^, we note that its usage leads only to 
qualitative description of the experimentally observed asymmetry of the scaling function. 
As can be seen from Figs. [914TT] the asymmetry due to the higher tail of f{ip) for positive 
values of ip is not very different for the three cases considered, the RFG model, the case with 
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HO single-particle wave functions and that with natural orbitals. In addition, the obtained 
asymmetry is rather similar in the cases when nonrelativistic or relativistic optical potentials 
are used to account for the FSI. Here we would like to emphasize that considerably larger 
asymmetry is obtained in the RMF+FSI method (see Figs. [T^ and [T^ . and especially when 
experimental widths with the Lorentzian function are used in the calculations (see Fig. [T^ . 
So, in the case of the RMF+FSI method the values of f{ip) for positive values of ip describe 
the experimental data quite satisfactorily. 

To summarize, our extensive study with different nuclear models based on the indepen- 
dent particle shell model and even incorporating short-range nucleon correlations through 
the use of NOs within the Jastrow method, show different results for the momentum dis- 
tributions, particularly at high missing momenta where nucleon correlations give rise to 
the presence of an important tail. However, the scaling functions obtained from them do 
not present the strong asymmetry at positive ip as given by the analysis of data. We have 
proven that such a "strong" asymmetry in f{ip) does not emerge even when a broader energy- 
dependence is assumed in the spectral function, but it does emerge when FSI are taken into 
account. In particular, taking into account final-state interactions in the RMF approach 
produces, in addition to a strong tail at high-p in the nucleon momentum distribution, a 
scaling function with the right amount of asymmetry in accordance with data. 

Here we would like to make some concluding remarks on the results obtained in the 
present work. We established and considered the relationship between the longitudinal scal- 
ing function f{ip) and the spectral function S{p,S) accounting for the restrictive condition 
of the PWIA in the description of the reaction mechanism and, correspondingly, with the 
momentum distribution at the specific conditions for the kinematically allowed region (the 
excitation energy Sraax — ^ oo). We studied the ingredients of the spectral function, namely 
the single-particle momentum distributions ni{p) and the hole-state energy distributions. In 
several consequent steps (starting with the IPSM) with increasing complexity beyond the 
MFA we showed that the ni{p) have to be considered within models in which the nucleon 
correlations are taken into account. This leads to the existence of high- momentum com- 
ponents of the momentum distributions. The methodical study of the energy distribution 
of the hole states showed that distributions with broader energy dependence (due to the 
residual interaction) have to be considered, with single-particle widths which are close to 
the empirically observed ones. The energy distributions must have more complicated form 
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than the Gaussian one, e.g. the Lorentzian one. However, it was pointed out that all men- 
tioned conditions are clearly necessary, but are not sufficient for a successful description of 
the experimentally observed asymmetry of the scaling function fXip), namely they give a 
weaker asymmetry at positive values of ip. We pointed out that this characteristic feature of 
f{ip) can be reached only if the FSI (and other peculiarities of the electron scattering beyond 
the PWIA) are carefully taken into account. This was reached in the case of RMF+FSI, 
while other methods using complex optical potentials for FSI lead to results that have some 
discrepancies with the experimental longitudinal scaling function f{ip) (e.g., for ip < —0.9 
and tp > 1.2), and they must be used with caution. The obtained results open a new 
perspective on the extraction of information on momentum distribution from experimental 
inclusive (e, e') data. 
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